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Abstract 



m 

J> ■ We present a heuristic, semiphenomenological model of the anomalous tern- 
(N 

^ ■ perature (T) dependence of resistivity pxx recently observed experimentally 

■ in the quasi-one-dimensional (QID) organic conductors of the (TMTSF)2X 

00 : 

0^ . family in moderately strong magnetic fields. We suggest that a QID conduc- 



tor behaves like an insulator {dpxx/dT < 0), when its effective dimensionality 



. is one, and like a metal {dpxx/dT > 0), when its effective dimensionality is 

O ■ greater than one. Applying a magnetic field reduces the effective dimension- 



ality of the system and switches the temperature dependence of resistivity 
between the insulating and metallic laws depending on the magnitude and 
orientation of the magnetic field. We critically analyze whether various mi- 
croscopic models suggested in literature can produce such a behavior and 
find that none of the models is fully satisfactory. In particular, we perform 
detailed analytical and numerical calculations within the scenario of magnetic- 
field-induced spin-density-wave precursor effect suggested by Gor'kov and find 

that the theoretical results do not agree with the experimental observations. 
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Typeset using REVTgX 



1 



I. EXPERIMENTAL INTRODUCTION 



In recent experiments a very unusual temperature (T) dependence of resistivity 

was observed in quasi-one-dimensional (QID) organic conductors (TMTSF)2C104 (at the 
ambient pressure) and (TMTSF)2PF6 (at a pressure about 9 kbar) in moderately strong 
magnetic fields H of the order of 10 T at T ^ 50 K. Unexpectedly large magnetoresistance 
in these materials has already attracted attention in earlier measurements 

The (TMTSF)2X materials consist of one-dimensional (ID) conducting chains parallel to 
the crystal axis a (for general reviews of the (TMTSF)2X materials, see Refs. [T|,^-|T2|). The 
chains are weakly coupled in the two others directions b and c, the coupling in the c direc- 
tion being much weaker than in the b direction. In zero magnetic field or in the field parallel 
to the b axis, the resistivity of (TMTSF)2PF6 along the chains, paa, depends on tempera- 
ture approximately quadratically: paa ~ T'^ |lllJ^JT3[|, which is consistent with the standard 
Fermi-liquid theory [|I4[, provided resistivity is dominated by electron-electron scattering. 
When a magnetic field is applied along the c axis, paa does not change appreciably at high 
temperatures; however, below a certain magnetic-field-dependent temperature Tmin ~ 20 K, 
resistivity starts to grow with decreasing temperature: dpaa/dT < at T < Tmin lilBil- 
In other words, the behavior of the system changes from metallic, dpaa/dT > 0, to insu- 
lating, dpaa/dT < 0, when the temperature is lowered below Tmin. The temperature Tmin 
increases with the increase of magnetic field. Such a behavior is very surprising in view of 
the fact that no thermodynamic phase transition is observed in the system at this relatively 
high temperature Tmin ~ 20 K. (A phase transition into the magnetic-field-induced spin- 
density-wave (FISDW) state takes place at the transition temperature Tc ~ 2 K, which is 
an order of magnitude lower than Tmin-) As the temperature is lowered further below Tmin, 
Paa(T) continues to grow until another temperature scale Tmax ~ 8 K < Tmin is reached. 
Magnetoresistance is huge at T ~ Tmax^ Paa aX H = 7.8 T is about 10 times greater than 
Paa at -ff = 0]. At the lower temperatures T < Tmax, behavior of the system starts to 
depend crucially on the exact orientation of the magnetic field P,p!5[]. If the magnetic field 
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lies in a plane formed by the direction a of the chains and another integer crystallographic 
direction, such as c or c + b, paa recovers the metallic behavior dp^a/dT > at T < Tmax- 
For other, generic orientations of the magnetic field, paa retains the nonmetallic behavior, 
either continuing to grow with decreasing temperature: dpaa/dT < 0, or saturating at a high 
constant value. The temperature Tmax does not depend appreciably on the magnetic field. 
If at T < Tjiiax the magnetic field is rotated in the plane perpendicular to the direction a 
of the chains, paa exhibits sharp minima when the field is aligned with the integer planes 
described above (the so-called "magic angles" effect). This effect was discovered experimen- 
tally earlier |[I6HT8| following the theoretical suggestion by Lebed' [|^,|^. Finally, at low 
temperatures, the system may enter the magnetic-field-induced spin-density-wave (FISDW) 
phase at the transition temperature ~ 1-2 K, where paa increases sharply. Other diago- 
nal components of the resistivity tensor of (TMTSF)2PF6 seem to behave similarly to paa'- 
The angular dependence of pcc at T < Tmax and the temperature dependences of pcc 



22,23] and p^b [p3H25[ appear to be qualitatively similar to those of paa- However, in the 



very recent measurements it was found that, when the magnetic field is oriented along 
the magic direction c, pcc monotoniuosly decreases with temperature, whereas paa and pbb 
exhibit a minimum at T = T^^^. When the magnetic field is tilted away from c, p^c develops 
a minimum at a temperature much lower than Tmin for paa and pbb- 

In (TMTSF)2C104, a magnetic field also causes paa{T), as well as the NMR relaxation 
rate 1/Ti, to grow with decreasing temperature at T < T^in 0111) which is an indication 
of a nonmetallic behavior induced by the magnetic field. On the other hand, it was shown 



in Refs. pT| , |27| , p8| that the angular and temperature dependences of pcc in (TMTSF)2C104 
are quite different from those in (TMTSF)2PF6 and can be interpreted within the standard 
Fermi-liquid picture. The behavior of (TMTSF)2C104 may or may not |^ be complicated 



by doubling of the crystal period in the b direction occurring in this material at T = 24 K. 
To avoid complications, we will focus our theoretical study on (TMTSF)2PF6, which has a 
simple crystal structure. 

The behavior of resistivity in (TMTSF)2PF6 completely contradicts the conventional 



Fermi-liquid picture of a metal with an open Fermi surface. In this picture, applying a 
magnetic field perpendicular to the direction of the chains should produce no or very little 
magnetoresistance, should not alter the metallic temperature dependence of resistance, and 
should exhibit no magic angles effect [^. Thus, unconventional theoretical approaches are 



required. 

II. HEURISTIC THEORETICAL PICTURE 

We suggest that the following theoretical picture may qualitatively account for the un- 
usual behavior of (TMTSF)2PF6 H. 



It is well known theoretically (see, for example, Ref. [0) that the orbital effect of 
a c-axis magnetic field if on a Q ID conductor is characterized by the cyclotron energy 
Eh = ebHvp/c, where e is the electron charge, b is the distance between the chains in the b 
direction, is the Fermi velocity, and c is the speed of light. For the realistic values of the 
model parameters (see Sec. |Vllip , we estimate that Eh/H ^ 1.8 K/T. The magnitude of 
the cyclotron energy. Eh ~ 14 K at if = 7.8 T, is close to the temperature of the resistivity 
minimum at that magnetic field, Tmin ~ 20 K. Taking into account that the minimum of 
resistivity clearly has a magnetic origin (it does not exist without magnetic field), and Tmin 
grows with the increase of magnetic field, we suggest that the minimum of resistivity occurs 
when the temperature reaches the energy scale of the magnetic field; that is, Tmin ~ Eh- 

Now we need to identify the nature of the second energy scale in the problem, the 
temperature of the resistivity maximum, Tmax- At the temperatures T > Tmax, it appears 
that paa depends only on the magnetic field projection on the c* axis perpendicular to the 



a and b directions From this observation, we may conclude that at T > Tmax the 



system behaves effectively as a two-dimensional (2D) system; that is, the coupling between 
the chains in the c direction is not relevant. On the other hand, at T < T^ax the coupling 
along the c axis becomes important. This is manifested by the magic angles effect, which 
is an essentially three-dimensional (3D) phenomenon involving both the b and c axes. The 
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coupling between the chains along the c axis is characterized by the electron tunneling 



amplitude tc, whose magnitude is believed to be of the order of 10 K ||33|, which is close to 
Traax ~ 8 K. Thus, wc suggcst that the electron tunneling amplitude tc sets the temperature 
scale Tjnax of the resistivity maximum: T^ax ~ tc. This conjecture is supported by the 
experimental fact that Tmax (unlike Tmin) does not depend appreciably on the magnetic field 
IP^. We also need to mention that, according to Refs. [|^,^, the coupling between the 
chains in the b direction, tb, is much greater than the temperatures discussed in our paper: 
tb ~ 300 K. 

Taking into account these energy scales, we identify three qualitatively different regimes 
in the behavior of a QID system in a magnetic field: 

1) High temperatures: Eh ^ T^m < T < tb- In this region, the temperature is greater 
than both the magnetic energy and the electron tunneling amplitude tc along the c axis, 
but lower than the tunneling amplitude tb along the b axis. Thus, we may neglect both the 
magnetic field and the coupling between the chains along the c axis and treat the system as 
a normal 2D Fermi liquid without magnetic field. This results in the quadratic law p^a ~ T"^ 
and the metallic behavior of the resistivity dpaa/dT > 0. 

2) Intermediate temperatures: tc ~ T^ax < T < T^nin ~ Eh- In this region, the temper- 
ature is still greater than the coupling between the chains along the c axis, so the system 
remains 2D; however, the effect of the magnetic field becomes important. It is known that, 
in the presence of a magnetic field along the c axis, the motion of electrons along the b 
axis becomes quantized, and the dispersion law of electrons becomes one-dimensional (ID) 
p6| . The degeneracy of the electron spectrum in the b direction is a specific manifestation 
of the Landau degeneracy in a magnetic field in the case of a 2D system with a strong 
QID anisotropy. This phenomenon is called "one-dimensionalization" of a QID system by 



a magnetic field |TT[]. Even though the spectrum of electrons becomes ID and their wave 
functions become localized in the b direction, the wave functions still spread over many 
chains (if Eh ^ t^), which results in a considerable interaction between different chains 



| 3T| . So the system is not truly ID, because it does not consist of completely decoupled ID 



chains. Nevertheless, we may expect that, at least, some ID features would be present in this 
regime and, via a mechanism that need to be specified, would lead to an insulating transport 
behavior dpaa/dT < 0. In general, ID systems have stronger tendency toward insulating 
behavior than higher-dimensional systems, because various insulating mechanisms, such as 
renormalization of umklapp scattering, density-wave instabilities, and Anderson localiza- 
tion, are more effective in one dimension than in higher dimensions. So the conjecture that 
the insulating behavior is caused by the magnetic-field-enforced "one-dimensionalization" 
is plausible, but requires detailed studying of a specific mechanism We review possi- 
ble candidates for the mechanism in the next section and quantitatively analyze one of the 
mechanisms in rest of the paper. 

3) Low temperatures: T < T^ax ~ ic- In this region, the coupling between the (a, b) 
planes becomes important. The magnetic field pointing exactly along the c axis does not 
affect the electron motion along that axis. Thus, in addition to the magnetic-field-enforced 
ID dispersion law discussed in part 2), the system acquired an extra dispersion in the c 
direction and becomes effectively 2D, which results in a metallic, Fermi- liquid behavior 
dpaa/dT > 0. If the magnetic field does not point along the c axis, the component of the 
field perpendicular to the c axis suppresses the energy dispersion along that axis, so the 
system remains effectively ID and insulating: dpaa/dT < 0. If the direction of the field 
is close to the c axis, we expect resistivity to decrease with decreasing temperature in the 
range E^^^ < T < Tmax ~ and to start increasing again at T < E^^\ where E^^^ is 
the cyclotron energy of the electron motion along the c axis, which is proportional to the 
projection of the magnetic field perpendicular to the c axis. The same arguments apply not 
only to the c axis, but also to the c + b axis and other integer crystallographic directions 
mc + nh. However, because the electron tunneling amplitudes in these directions decrease 
rapidly with the increase of the integers m and n, the effect is clearly visible experimentally 
only for the c -|- b axis. 

In summary, we suggest that the unusual transport behavior of (TMTSF)2PFe results 
from the changes in the effective dimensionality of the system caused by the applied magnetic 
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field. The system is 2D at Eh < T < th and effectively ID at tc < ^ < Eh- At T < tc the 
system is effectively 2D for the magic orientations of the magnetic field and effectively ID 
for generic orientations. Whenever the system is 2D (or 3D), it is a normal Fermi liquid, and 
the temperature dependence of resistivity is metallic. Whenever the system is effectively 
ID, the temperature dependence of resistivity is insulating. The latter state of the system 
might be called the magnetic-field- induced Luttinger insulator (MFILI), by analogy with 
the term "Luttinger liquid" , which refers to the metallic state of a ID system . 

We do not have detailed mathematical calculations that can prove the heuristic picture 
outlined in this section. Nevertheless, we can predict some experimental effects based on this 
picture. In Refs. p8|j3^ , oscillations of pcc upon rotation of a magnetic field in the (a, c) plane 
were discovered by Danner et al. Following the theoretical suggestion of Ref. it was 



found that a small magnetic field along the b axis destroys the oscillations [Q. We predict 
that if a magnetic field is rotated in the magic plane from the b + c direction toward the a 
direction, the Danner oscillations should exist, even though the magnetic field has a finite 
b-component. The suggested geometry has an advantage over the geometry of experiment 
where the magnetic field had a fixed b-component, that the Danner oscillations would 



not be mixed up with the Lebed' oscillations occurring when the magnetic field is rotated 
in the (b, c) plane. This prediction is based on the idea that the Danner oscillations require 
that the electron motion in the third direction is not suppressed by the magnetic field, which 
happens only when the magnetic field belongs to a magic plane. We also predict that the 
Danner oscillations should disappear at T > Tmax? where the electron dispersion in the third 
direction is smeared out by temperature. A detailed study of magnetic oscillations upon 
rotation of the magnetic field around the c* axis with different tilts relative to the c* axis 
was recently performed by Lee and Naughton They found possible to interprete most, 
but not all, of the results within a conventional semiclassical theory of metals. 
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III. REVIEW OF THEORETICAL MODELS 



Within the heuristic framework presented in the previous Section, a theoretical study of 
the problem reduces to the following two parts: 

(a) How a magnetic field induces the negative temperature dependence of resistivity 
dpaa/dT < in a 2D metal with a strong QID anisotropy. (2D problem) 

(b) How the electron tunneling in the third direction does or does not suppress the effect 
found in part (a) depending on the orientation of the magnetic field in the (b, c) plane. (3D 
problem) 

From the above formulation, it is clear that the 3D problem (b) can be addressed only 
after the 2D problem (a) has been solved. In other words, we believe that the insulating 
temperature dependence of resistivity (problem (a)) and the drop of resistivity at the magic 
angles (problem (b)) have a common origin. 

However, until recently, theoretical and experimental efforts were focused on solving 
problem (b) without recognizing and addressing problem (a). Soon after the experimental 



discovery of the magic angles |Tq-[T8| , a number of theories tried to explain this effect semi- 
classically |41|^3|. The theories [^,0 found the magic angles effect in p^b and pcc, but not 



in paa- Microscopic analysis I^O] of the "hot spots" model |^ demonstrated that it cannot 
explain the huge magnetoresistance and the magic angles effect. In all of these models, re- 
sistivity was calculated by introducing a relaxation time r phenomenologically and studying 
semiclassical electron trajectories in the momentum space. These theories assumed that 
relaxation mechanism does not change dramatically as the magnetic field is rotated. In view 
of experiment 0], this assumption is completely wrong, because it is the switching between 
metallic and insulating states and the corresponding change in the relaxation mechanism 
that causes the magic angles effect. 

Another explanation of the magic angles was proposed on the basis of the Luttinger 
liquid concept in Ref. ||3^ (see also Ref. |0). This theory suggested that the magic angles 
effect reflects the change of the effective dimensionality of the system: The system is a 2D 
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Luttinger liquid at generic angles, and a 3D normal Fermi liquid at the magic angles. The 
prediction of the theory that even a small b-component of the magnetic field would destroy 
coherence of the interchain hopping in the c direction was confirmed experimentally in Ref. 



3^,^. However, this theory focuses only on problem (b), but does not address the issue 
of the temperature dependence of resistivity and does not explain how problem (a) may 
be solved. Moreover, in the actual calculations ^9], tc is treated as a perturbation to a 



ID system, and the much greater tunneling amplitude tb is effectively neglected. For this 
reason, the theory |^ actually studies the dimensionality crossover between 2D and ID, 
not 2D and 3D, and the very important c-component of the magnetic field does not appear 
in these calculations. 

The very first calculation of the angular dependence of paa was done by Lebed' and Bak 
20| before the experimental discovery of the magic angles. However, it predicted maxima. 



not minima, of resistance at the magic angles. This discrepancy was corrected in the sub- 



sequent work [0. In this theory, resistivity p^a oc 1/r is studied by calculating the rate of 
umklapp scattering 1/r in the lowest order of perturbation theory. This approach allows to 
study how resistivity depends on temperature and on the magnitude and orientation of the 



magnetic field. According to Ref. PS|, the temperature dependence of the scattering rate 



changes with the magnetic field orientation: 1/r oc at the magic angles and 1/r oc T at 
generic angles, because the electron dispersion law is 2D and ID in these cases, respectively. 
At low enough temperatures, the difference between the and T laws should result in 
sharp dips of resistance at the magic angles. However, the theory predicts that both tem- 
perature dependences are metallic {dpaa/dT > 0), whereas experimentally the temperature 
dependence is insulating {dpaa/dT < 0) for nonmagic angles at T < T^in and for all angles 
at ^max < T < Tmin IQ-Q. Thus, the theory |^ is not adequate either. 



In Ref. Lebed' and Bak considered also the scattering rate of electrons on impu- 
rities renormalized by the electron-electron interaction in the lowest order and noted that 
it grows with decreasing temperature, when a magnetic field is applied. It is easy to check 
analytically that this scattering rate is approximately constant without magnetic field and. 



if tc is neglected, grows as In^En/T) in a magnetic field at T < Eh- We calculated this 
diagram numerically in the 2D case and found that it increases by only about 20% in the 
relevant range of fields and temperatures, which is insufficient to explain the experiment. 
(Our numerical calculation of the umklapp electron-electron scattering diagram in the 2D 
case also confirmed that application of a magnetic field changes the temperature dependence 
of 1/r from to T at T < Eh, as discussed in the preceding paragraph.) 

After the magnetic-field-induced insulating temperature dependence of paa was discov- 
ered experimentally, theory started to address this problem specifically. Refs. suggested 
that resistance increases at T < Tmin, because the "one-dimensionalization" of (TMTSF)2X 
by a magnetic field induces formation of a pseudogap in the charge channel. It is well known 
that charge and spin excitations are independent in a ID system [|12[. In the presence of 
umklapp scattering and repulsive interactions, charge excitations may develop a pseudogap, 
whereas spin excitations may remain gapless. This results in insulating temperature depen- 
dence of resistivity coexisting with metallic behavior of spin susceptibility. This effect is 
observed experimentally in the sulfur-based compounds (TMTTF)2X which are more 
ID than the selenium-based compounds (TMTSF)2X (see Refs. P,P|,P^). Induction of a 



charge pseudogap by a magnetic field would explain insulating behavior of all components 
of the resistivity tensor in (TMTSF)2PF6 l2|-|5| and the NMR data in (TMTSF)2C104 B 
(although the transverse resistivity in (TMTSF)2C104 behaves differently p7|p8| ). 

However, in order to achieve quantitative agreement with the experiment, this theory 
assumes that t;, ~ 50 K p[ or 15-30 K [^, which is almost an order of magnitude smaller 



that the commonly accepted value tj, ~ 250 K deduced from the b-axis plasma edge |35| and 



the (a,b)-plane angular magnetic oscillation In this theory, the magnitude of the pseu- 



dogap is determined by the energy scale where the renormalization-group (RG) equations 



(also called the parquet equations for the forward, backward, and umklapp scatter- 
ing amplitudes of electron-electron interaction diverge. However, because the chains are 
strongly coupled, it is not correct to limit the RG equations to only those three amplitudes 
of interaction. It is necessary to include an infinite number of the amplitudes of interaction 
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between different chains, which maizes the RG equations integro-differential |^. For a QID 
conductor in a magnetic field, the integro-differential RG equations were derived and solved 
numerically without umklapp in Refs. pl|j49|j50|] . The temperature where a solution of the 



RG equations diverges was interpreted in the latter papers as the FISDW transition tem- 
perature Tc. On the other hand, in a purely ID case, the energy where a solution of the RG 
equations diverges is conventionally interpreted as a pseudogap energy, not as a transition 



temperature, because thermodynamic phase transitions are not possible in ID systems [pi 
Thus, the problem is whether the RG equations for a QID conductor in a magnetic field can 
simultaneously describe formation of a pseudogap at a high temperature and the FISDW 
transition at a much lower temperature. In a purely ID case this is possible, because the RG 
equations separate exactly into two independent sets of equations for the spin and charge 
channels 1^^. However, in a higher-dimensional case, all of the interaction amplitudes are 
coupled, and separation of the RG equations into independent channels does not seem fea- 
sible. It is not clear why a pseudogap energy and a FISDW transition temperature would 
differ by an order of magnitude in the RG approach, if both are produced by the same 
mechanism of "one-dimensionalization" and enhancement of the Peierls susceptibility by a 
magnetic field. 

Even after the RG equations are solved, the temperature dependence of resistivity still 
needs to be calculated. This issue was addressed in the theory by Gor'kov which 



suggested that resistivity is given by the same diagram that was studied in Refs. pO| , ^ , but 
with a renormalized, temperature-dependent umklapp scattering amplitude. In the stan- 
dard model of the FISDW transition [|5^, where tf, ^ Eh, and electron-electron interaction 



is repulsive, it is safe to neglect the superconducting channel contribution to the RG equa- 
tions [^]. In this case, the RG equations reduce to the conventional ladder- RPA equations. 
The umklapp scattering amplitude, renormalized according to these equations, grows when 
T Tc. Because paa is proportional to the square of the umklapp amplitude, the insulating 
regime dpaa/dT < might be achieved if the umklapp amplitude grows faster than the 
phase space factor T or decreases. One might expect, though, that resistivity would 
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start to grow only in a narrow vicinity of the FISDW transition, not at the temperatures 
an order of magnitude higher than T^. Whether this approach can quantitatively explain 
the experiment, particularly the role of a magnetic field, can be verified by numerical cal- 
culations. In the rest of the paper, we study this issue in detail. The approach of Refs. 
5^,|55[] is attractive, because it permits straightforward calculation of transport coefficients 



and connects naturally with the standard, successful model of FISDW on one hand and with 
the simpler, better-understood transport model of Refs. [^,^] on the other hand. 



IV. QID CONDUCTOR IN A MAGNETIC FIELD 

For the following theoretical description, we select the x, ?/, and z axes along the crystal 
directions a, b, and c, which are not orthogonal in the triclinic (TMTSF)2X crystals. Elec- 
trons can tunnel between different chains with the amplitudes ti, where 1 = (Ix^l'i) is a 2D 
integer vector that determines the transverse displacement of the electron, di = /ib + Z2C, in 
the process of tunneling. The Fermi surface of ID electron motion along the chains consists 
of two Fermi points characterized by the Fermi wave vectors ifci?. We label the electrons 
with the wave vectors close to i/c^? by the index a = ±. In the vicinity of the Fermi 
energy, the energy dispersion law of the longitudinal electron motion can be linearized as 
e = ^hvpk^, where h is the Planck constant, vp is the Fermi velocity, the energy e is counted 
from the Fermi energy, and the longitudinal wave vector is counted from ^kp for the 
two Fermi points. In the absence of magnetic field, the total, longitudinal and transverse, 
electron dispersion law is 

£Q,(k) = ahvpkx + ^ 2ti cos(k_Ldi + aifi), (1) 
1 

where k = {kx, ky, k^) = {kx, kj_) are the electron wave vectors along the x, y, and z axes. 
The phases ipi in the transverse dispersion law (|1]) of a triclinic crystal are determined by 



the amplitudes of tunneling to different molecules belonging to the same chain |56]. In this 



paper, we present analytical formulas for a general dispersion law (|I]) with any number of 
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the transverse tunneling amplitudes t\, but we perform numerical calculations only for the 
2D case with the two tunneling amplitudes: tf, — ti Q between the nearest and — t2,o 
between the next-nearest chains in the b direction. The values of the corresponding phases, 
(pb and (p'f^, are not known reliably. According to Yamaji = — 7r/2, and ipb varies from 

7° to 40° when temperature varies from 300 K to 1.7 K in (TMTSF)2PF6. In our numerical 
calculations, we assume (f'f^ = 2(fb for simplicity and consider several values of (fb between 
the two extremal values (pb = and (pb = 7r/4. 

Suppose that a magnetic field H is applied perpendicular to the chains. It can be 
introduced into the Hamiltonian of the system via the Peierls-Onsager substitution k — > 
k — x[H X erc]e/hC, where e^, is the unit vector along the x axis, e is the electron charge, and 
C is the speed of light. The eigenfunctions ipa,\i of noninteracting electrons in the magnetic 
field H are ^ 



'ipa,k{^,ny,n^) = ^1— , exp 



' 2t\ 

kxX + kyUyb + kzTizC + a V - — — sin(k_Ldi - dx + aipi) 



(2) 



where Gi = ■ [di x Hje/hC are the wave vectors proportional to the magnetic field, ny 
and riz are the integer coordinates of the chains in the y and z directions, L is the length of 
a chain, and Af is the total number of the chains. The eigenenergies of eigenfunctions @ 
are 

e = ahvpkj., (3) 

thus the electron dispersion law is effectively ID in the magnetic field. 

The single-particle Green function of noninteracting electrons in the magnetic field H 
was found in Ref. p6 |: 

r°° dkx exp[iA;^.(x — x')] 



Qc,{x,x\'k^,iujm) = / 

J — ( 



27r iujm — ahvpkx 

X exp ( - — |^[sin(k_Ldi — G\x + a(p\) — sin(k_Ldi — G\x' + ay^i)] ] , (4) 
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where Um = 27r(m + 1/2)T is the Matsubara frequency. The Green function is a product 
of two terms: The first term is the Green function of ID electrons, whereas the second, 
exponential term represents the transverse motion of the electrons. Only the second term 
contains the magnetic field via the wave vectors Gi, which introduce periodic dependences 
on X and make the Green function (^) not translationally invariant along the chains. In 2D 
case, we denote the magnetic wave vector Gi^o as simply G. The integral over kx in Eq. 
(I) can be easily taken, but the Lehmann representation is more convenient for analytic 
continuation from the Matsubara frequencies to the real frequencies. 

V. RENORMALIZATION OF THE UMKLAPP SCATTERING AMPLITUDE 

DUE TO SDW INSTABILITY 

The tendency of a QID system toward the SDW instability manifests itself in divergence 
of the density-wave susceptibility, which is shown diagrammatically in the lowest order in 
Fig. |I[ In this figure, the solid and dashed lines represent the Green functions of the + and 
— electrons, and Qm and q_L are the incoming energy and the transverse wave vector. In 
the Matsubara technique, the analytic expression for the bare susceptibility per one chain 
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IS 



be (i^k_|_ 

T X{x' — x,q±) r°° dw sm.[{x' — x)w /Tivp] 



/be rC_|_ 
-J—^G^{x,x',^^^ + ^_^,^Uj + ^Vlrn)G+{x',x,^<i^,^u:j 

(5) 



2'k(1ivfY sinh[27rT(x' — x)/Tivf\ J-oo w — iVt. 



where 



A(x, q_L) = j exp ^-i ^ f~Q^ sin(G'ia;/2) cos(k^di) cos(q^di/2 - (/^i) j . 



(6) 



In Eqs. d^) and (^) and elsewhere, the integration over the transverse wave vectors k_|_ 
goes over the Brillouin zone. Because of the averaging over k^, susceptibility (^) is a 
translationally invariant function of x' — x, unlike the Green function (^). This allows us 
to Fourier-transform Eq. @ over x' — x and to obtain xo(ct;'^^m) as a function of the 3D 
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wave vector q. As follows from Eq. Xo(a;, q±, «f^m) behaves as 1/x when x < vpI^-kT. 
This results in logarithmical divergence of xo(q)^^m), which we cut off at a small distance 



xq = 1/2'ykF [Q, where 7 is the Euler constant ||55|| . 

Because the electron conduction band in (TMTSF)2X is half-filled, the Fermi wave vector 
is commensurate with the crystal lattice wave vector along the chains: 4k p = 27r/a. This 
relation permits the umklapp scattering process, where two "+" electrons are transformed 
into two "-" electrons, and the change of the total electron wave vector, Akp, is absorbed 
into the lattice wave vector 2Tx/a. The amplitude of this process is conventionally denoted 
by 73 1^31 • '^^^ one- loop diagram of Fig. generates a ladder renormalization of the vertices 



of interaction between electrons as shown diagrammatically in Fig. ^. In this figure, the 
wavy lines represent the bare, unrenormalized vertices of forward ((72) and umklapp {g^) 
scattering, the circles represent the corresponding renormalized vertices 72 and 73, and the 
thin lines inside the circles indicate spin conservations along the electron lines. Appropriate 
for the SDW channel, the + and — electron lines in Fig. ^have opposite spins. The equations 
of Fig. 1^ are the same as in Ref. [^. The analytic expression of these equations is 



72(q, iVtm) =92 + 5'2Xo(q, «f^m)72(q, i^m) + 5'3Xo(-q, -i^mhsifl, i^m) , (7) 

73(q, in^) =93 + fi'2Xo(-q, -^^m)73(q, i^m) + 5'3Xo(q, if^m)72(q, i^m), (8) 
and the solution is 

72(q i^m) = ^2 - i9l - 9i)Xoi-(l, -i^m) 

^ ' [1 -5'2Xo(q,^fim)][l -5'2Xo(-q, -^^m)] -5'3Xo(q,^fim)Xo(-q, -^^m)' 

73(q,in„) = y- — %r-T] 2 — ? — 7r~\ — ? tTT- ^^0) 

Notice that whenever the umklapp vertex appears in the r.h.s. of Eqs. (0) and (|^), the signs 
of q and iQm in the argument of Xo are reversed. 

As temperature decreases, the bare susceptibility Xo(q5 0) grows until the denominator 
in Eqs. (^) and ( PDj ) vanishes at a certain temperature Tc, which is the FISDW transition 
temperature: 

[1 - ^72Xo(q, 0)] [1 - ^72Xo(-q, 0)] - ^?|xo(q, 0)xo(-q, 0) = 0. (11) 
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The interaction amplitudes 72,3 (q,0) diverge at the transition temperature. Since Eq. ([TlD 
is a quadratic form of XolqjO) and Xo(~c[)0), it vanishes at two different temperatures for 
a given value of q. Usually, it is assumed that only the higher temperature is physically 
significant, and the transition temperature is determined from Eq. (|1T]) by selecting the wave 
vector q that provides the maximal value for T^. 

By considering equations similar to Eq. ([Tl|), Lebed' [^] predicted that the umklapp 



splitting of the FISDW instability would result in oscillations of Tc vs H, but the effect was 
not observed experimentally. On the other hand, the experimental T-H phase diagram of 
(TMTSF)2PF6 can be well reproduced while neglecting umklapp |^4[. In order not to spoil 
the phase diagram by the umklapp splitting and to avoid unnecessary complications, we 
assume that is sufficiently small and neglect it in Eq. ( p!TD and in the denominator of Eq. 
(P^) |58[| . Using the conventional band-structure parameters of (TMTSF)2PF6: ta = 2900 K 
(Ep = V^ta), tfe = ta/10, and vp = 2 x 10^ m/s, we find that the choice tj, = 20 K and ^2 = 
g2/2nhvF = 0.2288 produces a T-H phase diagram close to the one observed experimentally 



in Ref. |^9[ at 12 kbar. The phase diagram does not depend on (pb and is shown in Fig. |^. 
In this figure, various symbols indicate the integer number N of the quantized longitudinal 
wave vector of FISDW: qx = NG . Note that wide spacing in H between the points of 



our calculations would not allow us to observe the oscillations of Tc vs H , even if we took 
into account the umklapp splitting of the FISDW instability. If we set t[ = 0, the electron 
spectrum ([l|) acquires perfect nesting: e+{kx, ky) = —e-{kx, ky — {n — 2Lpb)/b), and the SDW 
transition temperature To = (/;,fi7'/7ra;o) exp(— 1/^2) = 14.7 K becomes independent of the 
magnetic field. 

While the Matsubara representation of 72,3 is useful for determining the T-H phase 
diagram, we need the scattering vertices at the real frequencies to calculate resistivity. An- 
alytically continuing Eq. from the Matsubara frequencies iQm to the real energies e, we 
find the following expression for the bare susceptibility: 

T f°° Re{exp{-iqxx)\{x,q±)}exp{ixe/hvF) 
Xo q,£ =77 — ^ dx . , rj. 1^ — ^ • (12) 

16 



Substituting Eq. (0) into Eqs. (|) and (|TOD, we find the scattering vertices 72,3(q, ^) at the 
real energies e. 



VI. UMKLAPP RESISTIVITY OF A Q ID METAL 



Using the variational principle for the Boltzmann equation we find the following 
expression for the resistivity along the chains due to electron-electron umklapp scattering: 



P^^=^-Wur ^ /o(^l)/o(^2)[l-/o(^3)][l-/o(^4)]<%'^ (13) 



where Ei = e(ki) are the energies of electrons in the eigenstates |kj) (0), fo{e) is the Fermi 
distribution function, and W^^'^^ is the scattering rate of two electrons from the states |ki) 
and |k2) into the states jka) and |k4): 



27J- 

= -J- I (ki, k2|73|k3, k4) P S{e, + 62-6^- e,). (14) 



Here 5{e) is the Dirac delta function, and the matrix element is 

/oo 
dxi dx2 V'+,ki(a;i, n^^^) iIj+^]^^{x2, nf^ ,nf^) 

yi) u; yz) 

xC,k3(3^2, nf,nf^) Ck4(^i. n^y\n^P) l^i^i - ^2, - nf ,n^P - , £2 - e^). (15) 

The vertex of interaction 73 is written in the mixed, coordinate and energy, representation, 
describing the umklapp scattering of the two "+" electrons with the coordinates xi and X2 
located on the chains {rif\ri{^) and (ng^"*, n2^'') from the states |ki) and |k2) into the two 
"— " states |k3) and |k4) with the same coordinates. 

Substituting Eq. (g) into Eq. (|1^) and changing the variables of integration k^^ — * ei via 
Eq. (^, we find the following expression for the resistivity 

2{-nnf L%c deide2 decide i w / //mm fi m w^a £4 (^a^ 

P^^ = g2y^ (2TxhvFY ~ /o(^3)][l - k{ei)] WI^J^, (16) 

where 
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£4 



E 



k3,k4 

£l,£2 /-^ ' ' ki,k2 



wt 



_ hcJ^ 
X exp 



dx' dx'' 



d^q 



(2vr)2 

i{qx + -j^ — I X 



|73(q, £2 - e^)\^5{ei + £2 - £^3 - £^4) 



A(x',-q_L)exp 



. , , ^2 + £^3 \ // 
I ( -qx + I X 



A(x",qx). 



(17) 



Using the 5-function from Eq. (0), we take the integral over £4 in Eq. ([16D . Then, changing 
the integration variables ei, £2, and £3 to e = £2 — £^3, £^'1 = £^1 + (£^2 — £^3)/2, and e'2 = £2 + £3 
and taking the integrals over e\ and e'g, we obtain the final expression for resistivity: 



where 



73 (q,^) 



32e2T J- 



2T 



/CO r 
de / rf^qF(q,e)F(-q,e) |73(q,£) 
-00 J 



(18) 



hvp sinh(e/2T) 

73 (q, 



dx 



exp{—iqxx)\{x, q±) sin^xe/hvp) 
smh.{27iTx/hvp) 



(19) 
(20) 



are dimensionless functions. The renormalized umklapp vertex 73(q, e), given by Eqs. (p!oD 
and (0), should be substituted into Eqs. (|18D and (pOl) . If the renormalization of 73 is 
neglected (73 = gs), the integral over e in Eq. ([T8| ) can be taken analytically, and the result 
agrees with the expressions obtained in Refs. pO,H5|. 



VII. TEMPERATURE DEPENDENCE OF RESISTIVITY IN ZERO MAGNETIC 

FIELD 

Having obtained the general expression (|TBp for the longitudinal resistivity, let us examine 
the limit of zero magnetic field first. In this case, F(q, e) ([19|) becomes 

T?f \ be r (i^k_L 

"''(^-J (^) + cosh (^"'^■'- + S' - '^"^"^^■J')) • ' ' 

The function F(q, e) (|2l|), confines integration in Eq. (|18D to the energy interval 

\e\ < T (22) 



and to the wave-vectors region defined by the two inequahties: 

Ihvpqx + Xl4^iCOs(q_Ldi/2 - (fi) cos(k_Ldi)| < T, (23) 
1 

Ihvpq, - E 4^1 cos(qxdi/2 + (^i) cos(k'^di) | < T. (24) 
1 

It follows from Eqs. (pSD and (^) that the integration over the transverse wave vectors is 
restricted by the inequality 

|S(q^,k^,kl)|<T, (25) 

where 

S(q±, k_L, k^) = 4ti[cos(q_Ldi/2 - (fi) cos(k_Ldi) + cos(q_Ldi/2 + (fi) cos(k'_Ldi)]. (26) 
1 

The integrals over e and in Eq. (0) with the function -F(q, e) (|2l|) can be taken 

analytically, provided we neglect dependence of 73(5'x, q±, £^) on q^^ and e, that is, replace 

Isiqx, q±, £) by 73(^1, q±, £^), where qx and e are some characteristic values of qx and e from 

the intervals of integration (|2|), (|2|), and (||): 

7r^(&c)^T f d^^d'k^d'c^^ ( S(q^,k^,kY)/r W 
^-"Te^y i sinh[S(q„k,,kl)/4T] j l^3(g.,q.,.)| . (27) 

This approximation should be valid provided the peak in 73(^1., q±, £^) at e = and q equal 
to the nesting vector is wider than temperature, which is the case for temperatures not very 
close to the transition temperature Tc. Eq. (|2^) with S (^) coincides with the expression 
obtained by Gor'kov for ip\ = 0, except that we find an additional factor 1/2 in the 



argument of sinh, as in Ref. [ |6T[] . 

If renormalization of 73 is neglected, and 73 is replaced by a constant gs, then integral 
(P7|) gives the volume of the wave- vectors space restricted by inequality (^). In a general 
2D or 3D case, this volume is proportional to T, thus resistance (^) is proportional to T^: 

^(2D,3D)^~2y2^ (28) 



which is the standard result of the Fermi-liquid theory |T^. In ID case |^, where S = 0, 
Eq. ( p7|) reproduces the result of Ref. : 
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pi^°^-^l7rfe,^-)P- (29) 



The ID resistance (pQ]) is proportional to temperature T multiplied by the square of the 
renormalized umklapp amplitude 73, which may also depend on temperature. 

Now let us consider the simplest 2D case where only one tunneling amplitude ti, is kept 
in Eq. (^). If y^j, = (similar equations hold also for 9?;, = 7r/2), then Eq. can be 
factorized [B2|,B3: 



S(g^, ky, k'y) ^^^^ = 8hcos{qyb/2) cos[{ky + k'y)b/2] cos[{ky - k'y)b/2]. (30) 
If renormalization of 73 is neglected (73 = g^), then Eq. (|27|) with E from Eq. (^) gives: 



.(2D) 
rxx 



^~giT'\n\U/T), (31) 



which has an extra logarithmic factor compared to Eq. ( pH] ) ||61|| . Eqs. ( p8D and (^) are in 
agreement with the results of Ref. where the so-called "hot spots" in the distribution of 
the umklapp scattering time over the Fermi surface of a QID metal were studied. The "hot 
spots" are the points where the scattering rate is strongly enhanced compared to the rest 



of the Fermi surface, typically by the factor ln{tb/T) and occasionally by the factor ^tb/T. 
Positions of the hot spots are determined by the saddle points of the function H(gj^, ky, k'y) 
| 29| . When Lpb 7^ 0,7r/2, only isolated hot spots exist on the Fermi surface. Because they 



occupy a small phase space, the hot spots do not contribute significantly to resistivity in 
this general case, and Eq. (|28|) holds. However, in the special case iph = 0,n/2, the entire 
Fermi surface becomes "hot" p9[, and resistivity acquires the logarithmic factor of Eq. (pll). 

In the special 2D case with only one tunneling amplitude tb and (pb = 0, the electron 
dispersion (|1]) has a perfect nesting at the wave vector = and Qy = ir/b. The zero-field 
susceptibility Xo(ct)^)) given by Eqs. ( [T^ ) and (|^), diverges logarithmically at e = = 
and qy = n/b. Since we neglect in the denominator of Eq. (pUj), the renormalized vertex 
73(0,7r/6, 0) becomes 

73(0,7r/&,0)oc (32) 
gi In (T/To) 
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where To = hvp / {ttxq) exp{~l /g2) is the SDW transition temperature. At the same time, 
condition (^) with S given by Eq. ( |30D restricts integration in Eq. ( pTj) to the vicinity of 
either qy = n/b or ky + k'y = n/b or ky — k'y = n/b. The first of these conditions is satisfied 
at the same wave vector qy = n/b where 73 (^) diverges. Assuming that the integral in Eq. 
(P7|) is dominated by the vicinity of qy = n/b, we find: 



.(2D) 
rxx 



vt=o 5,4 In (T/To) 

Eq. (|33|) is analogous to the expression obtained by Gor'kov [^,^, but differs in the powers 



of the logarithms. The factor in Eq. (|33|) tends to diminish resistivity with decreasing 
temperature, which is characteristic for a metal. On the other hand, the logarithmic factors 
in Eq. (0), both in the numerator and denominator, tend to increase resistance, which is 
characteristic for an insulator. Which of these two competing tendencies wins can be found 
numerically. 

In Fig. H, we show temperature dependences of resistivity PxxiT) calculated via Eqs. (|T8|) 
and (|21| ) in 2D case at zero magnetic field. Because we neglected in the denominator of Eq. 
(p!0|), pxx oc g1 exactly. In Figs. ||, |, and|^, we plot the ratio Pxx/g1 where the dimensionless 
umklapp scattering amplitude ^3 = g^l2T:Tivp cancels out. As discussed in Appendix, the 
value of ^3 can be recovered by comparing these figures with the experimental data. The 
top solid curve in Fig. ^ shows pxx{T) calculated with only one tunneling amplitude and 
ifh = 0. We observe that PxxiT) has a positive, metallic slope dpxx{T) / dT > at high 
temperatures, whereas at lower temperatures the slope is negative. It is surprising that the 
negative slope starts at the temperature about 60 K, which is several times higher than the 
SDW transition temperature Tq = 14.7 K. This indicates that the logarithmic factors in Eq. 
(p3|) overcome the factor at relatively high temperatures. 

Now let us discuss a more general 2D model with only one tunneling amplitude t^, 
but with ifb 7^ 0. In this case, 73 diverges at the nesting vector = and qy = {ji — 
2ipij)/b. The umklapp vertex 73(0, (tt — 2(y9b)/6, 0) oc 5'3/ln(T/To) diverges at exactly the 
same transition temperature Tq as in the case = 0, but less strongly than in Eq. (^2|), 
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because Xo(0, qy, 0) and Xo(0, —qy, 0) in Eq. (|TOD do not diverge at the same wave vector qy. 
At ipb 7^ 0, condition (|25| ) with H given by Eq. ( |26D is not satisfied at the nesting vector 
% = ~ 2(y9fe)/6 independently of ky and ky, which further reduces pxxiT) compared to 
Eq. (|33| ) for (^i, = 0. Temperature dependences PxxiT) are shown in Fig. || by sohd curves 
for different values of y^f,. All curves have the same values of tf, = 290 K, g2 = 0.2288, and 
To = 14.7 K as the top solid curve with ipb = 0. While all solid curves diverge at the same 
transition temperature Tq, the region of the negative slope in pxx{T) shrinks rapidly with 
increasing ipi, and becomes much smaller than Tq at (fb ^ 7r/27. Thus, a non-zero phase ifb 
strongly suppresses the precursor effect in resistivity. For the curves with ~ 7r/27, the 
behavior of Pxx{T) qualitatively follows the 2D quadratic law (^) at the low temperatures 
T < 2t(,/7r = 185 K and the ID law (^) at the higher temperatures T > 21^1 t^- The slope of 
Pxx{T) at T ^ 2tfc/7r is rather small, presumably because the ID logarithmic renormalization 
of 73 partially compensates the linear temperature factor in Eq. (p9|). 

The dotted curves in Fig. ^ represent Pxx{T) for the 2D model where the tunneling 
amplitude to the next-nearest-neighboring chains, = 20 K, is introduced in addition to 
tfe. Since t'^ eliminates nesting in the dispersion law (|I]), the system does not have SDW 
instability provided t[ > Tq, so 73 does not diverge. The slope of the curve with ipb = 
changes from negative to positive at T ^ '^t'b/iT, which creates a maximum in the Pxx{T) 
curve. Thus, cuts off the precursor effect in pxx{T) at T ~ as discussed in Refs. 



52,53 



In conclusion, we have confirmed the suggestion by Gor'kov that, in zero magnetic 

field, the renormalization of umklapp amplitude due to proximity to a SDW transition can 
produce a negative slope in Pxx{T) at temperatures much higher than the SDW transition 
temperature Tq. This precursor effect is suppressed by a non-zero phase in the electron 
dispersion law, which does not influence the SDW transition temperature, but shrinks the 
temperature region of the negative slope. The second tunneling amplitude tj, suppresses the 
SDW transition temperature and cuts off the negative slope at T ^ tyn. In (TMTSF)2X, 
negative slope of pxx{T) is not observed in zero magnetic fleld and appears only when 
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a magnetic field is applied. According to the scenario suggested by Gor'kov ||5^j53|] , the 
negative slope is eliminated by a non-zero t'j, (and, possibly, by a non-zero y?;,, as shown 
above), but it is restored by a magnetic field. We examine feasibility of this scenario in the 
next section. 



VIII. TEMPERATURE DEPENDENCE OF RESISTIVITY IN A MAGNETIC 

FIELD 

In 2D case, the effect of a magnetic field on a QID electron system is characterized 
by the wave vector G = ebH/hc and the cyclotron energy Eh = hvpG = ebHvp/c ||64|| . 
Using the parameters vp = 2 x 10^ m/s and b = 7.7 A, we find that Eh/H ^ 1.8 K/T. 
Magnetic field enters Eqs. (|1^) and (|19D for resistivity only through the function X{x,qy) 
(^). In 2D case, magnetic field makes A a periodic function of x with the period 27r/G: 
X{x, Qy) = \{x + 2-n /G, qy), so it can be expanded into a Fourier series: 

oo 

\{x,qy)= ^n(g,)e^"^^ (34) 

ra=— oo 

with some coefficients An- Substituting Eq. (^^ into Eq. (|l^), we find: 

£ . '''•'!°L ,„ _ „nu ■ (35) 



n=-oo cosh (^) + cosh ^M^i^—Jl^'j 



As follows from Eq. (|35|) , the integration in Eq. ([TsD is concentrated in the energy interval 
\e\ ^ T and, at T <^ Eh, in the vicinity of the integer wave vectors ~ nG. Assuming that 
the characteristic width in q^ and e of the function 73 (q, e) is greater than temperature, one 
could replace 73(q, e) by 'ys{nG,qy,0) in the integral ([T8|) . In this case, the integrals over e 
and qx can be taken analytically: 



Pi'.''^ - ^5^^ / d% £ \UnG,qyM 



00 

X E An,{qy)An^{-qy) [ . ' '"^ "^Tf^/^^i 1 • (36) 

ni,n1^-oo Vsinh[(ni + n2-n)Ej^/4T]y 



(m + na - n)EH/4:T 



2 
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At temperatures much lower than the cyclotron energy, T ^ Eh, only the term with 
Til + n2 = n contributes significantly to Eq. (R3): 



„ oo oo 

J dqy J2 \^3{nG,qy,0)\^ ^ An,{qy)An-mi-qy)- (37) 



Eq. ( P7D is similar to the ID formula (pUj), except for the additional integration over qy and 
summation over n. This is a consequence of one-dimensionalization of the electron spectrum 
(I) by the magnetic field. 

If renormalization of 73 is neglected (73 = gs), then Eq. (|37|) produces a linear tem- 
perature dependence for resistivity in agreement with Lebed's results p0| , |45[| for T ^ Eh- 
When renormalization of 73 is taken into account, 73(q, 0) diverges at a certain wave vector 
q^ = NG and qy = Qy as T Tc. Assuming that only the term with q^ = NO dominates 
the sum ([37|), we find: 



XX 2 



dqy\UNG,qy,0)\^ ^ A4qy)AN_^{-qy). (38) 



It is important to emphasize that the function -F(q, e) (p5D, while restricting integration over 
e and qx, does not restricted integration over qy significantly (unlike in zero magnetic field, 
Sec. |Vli| ), because the Fourier coefficients An{qy) are nonsingular, temperature-independent 
functions of qy. Thus, the integrals over qy in Eqs. (|36D -(p8l) are not restricted to the vicinity 
of the FISDW vector Qy. This does not allow us to replace •js^NG, qy, 0) by •y^^NG, Qy, 0) 
and take the latter outside of the integral, like in Eq. (|33D . The integration over qy reduces 
the divergence of Eq. ( ^Sf ) at T ^ Tc and makes the resistivity precursor effect of SDW 
in a magnetic field weaker than without the field. This happens because the phase-space 



restrictions discussed in Sec. |VII| are consequences of the 2D nature of the electron dispersion 
law ([1|) at zero magnetic field, whereas a nonzero field makes the electron spectrum (^ one- 
dimensional. 

In Fig. H, we show temperature dependences of resistivity calculated via Eqs. (|T8|) and 
(P5|) for t'f^ = 20 K, the phases iph = 0, vr/27, and 7r/4, and the magnetic fields H = 0, 5, 15, 
and 25 T. In the case of = (the top panel), only a very strong magnetic field H = 25 
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T restores the negative slope of resistivity, dp^x/dT < 0. This result cannot be applied to 
explaining experiments because for (pi, = the slope of resistivity is already negative 

at 15 < T < 60 K in zero magnetic field, which does not agree with the experiment. On 
the other hand, once we increase ipb to make the zero-field slope positive, the effect of the 
magnetic field becomes very weak outside of a narrow vicinity of for both a very small 
phase (fb = vr/27 (the middle panel of Fig. |^) and a rather big phase 7r/4 (the bottom panel 
of Fig. H). If we keep (fb = and increase t'f, to the values 30 K or 40 K, the zero- field slope 
becomes positive, but only an enormous magnetic field of 50 T makes the slope negative 
(see Fig. [^. One may conclude that in the case (pb = the negative slope of resistance 
occupies a substantial range of magnetic fields only at very strong fields such that > t'^: 
if = 25 T for = 20 K (the top panel of Fig. |) and ii = 50 T for = 30 and 40 K (Fig. 1). 
This scenario is hard to reconcile with the experiment, because, according to the standard 



theory of FISDW 0], the condition Eh > corresponds to the magnetic fields where the 
last, = phase transition in the FISDW cascade takes place, whereas the negative slope 
in resistance occurs at much lower magnetic fields. 

We conclude that the FISDW precursor scenario ||5^ , |53| cannot explain the anomalous 
temperature dependence of resistivity in QID conductors in a magnetic field observed in 
experiments • The insulating temperature dependence of resistivity, dp^x/dT < 0, does 
exist in this scenario at T ^ for a certain choice of the transverse dispersion law of 
electrons, however the effect is present even without magnetic field in this case. If the 
transverse dispersion law is modified to suppress the insulating behavior in zero field, then 
applying a magnetic field produces a negative slope in the temperature dependence of Pxx(T) 
only either in a narrow vicinity of Tc or at very strong fields such that Eh > t'^,. 
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IX. TEMPERATURE DEPENDENCE OF NUCLEAR MAGNETIC 

RELAXATION 



In Ref. [0, temperature dependence of the NMR relaxation rate in (TMTSF)2C104 was 
measured and was discussed as an evidence for a charge pseudogap formation. In this 
section, we study the effect of a magnetic field on the NMR relaxation rate within the 
FISDW precursor scenario. 



According to the theory of the NMR relaxation via the electron spin fluctuations ||66 
(see also Ref. [^), the NMR relaxation rate 1/Ti is proportional to the imaginary part of 
the electron spin susceptibility: 

1^^ lmx(*i^_ (39) 

J-1 q 



where Un is the nuclear Larmor frequency. The coefficient of proportionality in Eq. (|39D 
depends on the hyperfine interaction parameters. In QID compounds, the sum over the 
wave vectors q in Eq. (|39|) can be separated into the contribution from the uniform spin 
susceptibility at close to = and the contribution from the antiferromagnetic spin 
fluctuations at close to qx = ^kp- 

Tf 1 = Tf 1 (g, ^ 0) + T,-\qx ^2kF). (40) 

In Refs. HJ67|], the temperature dependences of each term in Eq. were measured sepa- 
rately in (TMTSF)2C104 in the magnetic field 15 T along the c* axis. It was found that the 
antiferromagnetic term T{~^{qx ~ starts to grow below the same temperature Tmin that 
separates the metallic and insulating temperature dependences of resistivity {dp^x/dT > 
and dpxx/dT < 0). It was claimed that the growth of T^^{qx ~ 2kp) is a manifestation of 
the opening of a charge pseudogap in the spectrum of electron excitations. 

In the ladder approximation, the renormalized spin susceptibility x(ct; t^) is given by the 
Feynman diagrams similar to those shown in Figs. |] and |^: 

xM = r, Xo(q,^)[l-«.Xo(-q,-^)] 

[1 -^2Xo(q,w)J[i -^2Xo(-q, -^>\ -^3Xo(q,w)xo(-q, 
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Neglecting the umklapp splitting of the FISDW instability, i. e. setting (73 = in the denom- 
inator of (1^), and taking the zero-frequency limit in Eq. (0), because the nuclear Larmor 
frequency is small compared to all other energies, we find: 

Ti^iqx ~ 2kF) given by Eq. (|4^ ) does not depend on the phase ipb- The temperature 
dependence of T^^{qx ~ 2kF) calculated via Eq. (|4^) with = 20 K is shown in Fig. 
0. We see that Ti^(qx ~ 2^^?) behaves in a magnetic field in the same way as resistivity 
does, i. e. it grows and deviates from the zero-field curve only in a narrow vicinity of the 
critical temperature. Thus, we conclude that the FISDW precursor scenario does not agree 
with the experimental behavior of the NMR relaxation rate in (TMTSF)2C104 in a strong 



magnetic field. The results of our calculations are very similar to those of Ref. |£8[ for a 
SDW transition without magnetic field [p5|]. 



X. CONCLUSIONS 

In this paper we presented a heuristic, semiphenomenological explanation of the anoma- 
lous temperature dependence of resistivity of QID conductors in a magnetic field observed 
in experiments reviewed in Sec. |. According to this picture (Sec. |I|), a QID conductor 
behaves like an insulator [dp^x/dT < 0), when its effective dimensionality is one, and like 
a metal {dpxx/dT > 0), when its effective dimensionality is greater than one. Applying a 
magnetic field reduces the effective dimensionality of the system and switches the temper- 
ature dependence of resistivity between the insulating and metallic laws depending on the 
magnitude and orientation of the magnetic field. Using this picture, we predicted that the 
Danner oscillations of pcc may be observed when a magnetic field is rotated in the magic 
plane from the b + c direction toward the a direction. We critically analyzed whether various 
microscopic models suggested in literature can justify our heuristic picture and found that 
none of the models is fully satisfactory (Sec. [m| ). We studied the FISDW precursor scenario 
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suggested by Gor'kov p2| , |53| in detail both analytically and numerically (Sec. ^Vl-p^) and 



found that it does not agree with the experimental observations. In the rest of this section, 
we speculate about possible alternative approaches to solving the problem. 

1. Within the FISDW precursor scenario, we studied only the umklapp scattering rate 
1/r. In the language of Feynman diagrams, 1/r is related to the imaginary part of the 
electron self-energy S. Via the Kramers-Kronig relations for E, any precursor effects in 
ImE should also affect ReE. Because ReS is related to the electron density of states, 
a pseudogap may open in the single-electron spectral density as a precursor of FISDW. 
Potentially, the pseudogap may affect transport properties of the system as strongly as the 
renormalization of scattering rate studied in this paper. ReS also determines the residue Z 
of the single-electron Green function, which is finite for a Fermi liquid and vanishes for a 
Luttinger liquid. ReS can be straightforwardly calculated using the methods of this paper. 
On the other hand, relationship between the single-electron and transport properties is not 
straightforward. For example, the residue Z may cancel out from resistivity due to the Ward 
identities [^. Reliable calculations of resistivity require taking into account corrections to 
the vertex of interaction between electrons and electromagnetic field, which is a difficult 
problem for inelastic electron-electron interaction [^. 

2. The ladder approximation utilized in this paper does not take into account contri- 
butions from superconducting channel, which, on one hand, plays a very important role 
in one dimension and, on the other hand, is sensitive to a magnetic field. The magnetic 
field may induce an insulating behavior by suppressing superconducting fiuctuations, which 
compensate insulating fiuctuations at zero magnetic field. Including both density-wave and 
superconducting channels naturally leads to the parquet equations. If the parquet equations 
in a magnetic field are generalized to include the umklapp amplitude, one may try 
to find out whether they can be decomposed into two disconnected sets of equation anal- 
ogously to the spin-charge separation in ID case. This is a nontrivial possibility, because 
the number of coupled equations is infinite, and the equations are nonlinear. Nevertheless, 
if the equations do decouple, the two sets of equations would diverge at two different tem- 
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peratures, one of which could be identified with Tmin and another with Tc. An additional 
technical problem is that the equations of Refs. pl|,^,|50| are applicable only at T ^ Eh, 



whereas Tmin ~ Eh, so the parquet equations need to be derived in the difficult range 
T > Eh- The parquet approach may require to assume a rather small effective value for th 
|7T|,|72|. 

3. Electron-electron scattering contributes to electrical resistivity only via the umklapp 
processes, which do not conserve the total electron momentum |60|. For commensurate 
systems, such as (TMTSF)2X, the umklapp scattering process that changes the total electron 
momentum by Akp = I'nja^ where a is the lattice spacing in the a direction, is usually 
considered. However, in the presence of a magnetic field, the total electron momentum may 
also change by a multiple of G, the wave vector of the magnetic field. Since this kind of 
umklapp exists only in a magnetic field, it may be a natural source of magnetoresistance 
in a QID system. It is reasonable to expect that this mechanism works effectively only at 
T < Eh- This could generate Tmin ~ Eh without invoking pseudogaps originating from 
renormalization, whose values tend to be close to Tc. On the other hand, this mechanism 
would not explain the insulating behavior of the transverse resistivities Pyy and pzz and 
the NMR relaxation rate. However, the recent experiment [E^, where metallic temperature 



dependence was found for p^z and insulating for p^x and pyy, may be in favor of this or 
another kinetic mechanism and against a charge-gap scenario. 

Undoubtedly, the anomalous temperature dependence of resistivity in QID conductors 
in a magnetic field poses a tough theoretical challenge and solving this difficult puzzle would 
greatly enrich the condensed-matter physics. 
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APPENDIX 



The absolute values of p^x in (TMTSF)2PF6 can be found in Figs. 1 and 2 of Ref. 0, 
Fig. 20 of Ref. [|], Fig. 4 of Ref. |T3[, and Figs. 2 and 3 of Ref. From these data, as well 



as from Ref. ||73|, we find that pxx ~ 18 pVt cm at T = 20 K. Comparing this experimental 
value of resistivity with the theoretical curves in Fig. ^, we find = 0.17 at = or 
^3 = 0.77 at = 7r/4. On the other hand, taking the value p^x ~ 1 cm at T = 300 K 
PJ^, we would find ^3 ~ 1. 

In Ref. [^], the experimental temperature dependence of pxx in (TMTSF)2PF6 was fit 



neglecting renormalization of 73, and the value g^^^^^ = 0.21 was found. However, because 
our Eq. (0) differs from Eq. (3) of Ref. by a factor {2tt)^, and a factor 0.5 is missing 
in Eq. (7) of Ref. |61[|, the value g^^'^^^ = 0.21 corresponds to g^ = 2Tx^/'^g^^^^'^ = 2.34 in 



our notation. The ratio of this value and the values quoted in the preceding paragraph, 
2.34/0.77 = 3.3, represents the effect of renormalization of 73 by the SDW diagrams. 
According to the Drude formula, resistivity along the chains is 

nbch 4:71 , . 

Pxx = —2 — = y^^) 

where r is the relaxation rate, and tu^ = Se^vp/bch is the plasma frequency. Using the values 
of vp, b and c, we find ujp = 970 cm~^, which is close to the value Up = 1.1 x 10^ cm~^ from 
Ref. Comparing Eq. ( ^31) with the quoted above experimental values of pxx, we find 



1/r ^ 40 K at T = 20 K and 1/r ^ 2200 K at T = 300 K. One may notice that both of 
these values for 1/r are greater than the corresponding temperatures: 1/r > T. For that 
reason, it may be more appropriate to use 1/r, rather than T, as an infrared cutoff of the 
renormalization. 
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FIGURES 

FIG. 1. Feynman diagram for the bare density-wave susceptibility Xoi^' ~ x,ci±, i^m)- The 
soUd and dashed hnes represent the Green functions of the + and — electrons. 

FIG. 2. Feynman diagrams for the renormalized vertices of forward (72) and umklapp (73) 
scattering in the ladder approximation. The wavy lines represent the bare, unrenormalized ver- 
tices of forward ((72) and umklapp (5(3) scattering, whereas the circles represent the corresponding 
renormalized vertices. The + and — electron lines (the solid and dashed lines) are implied to have 
opposite spins. Spin is conserved along the thin lines inside the circles. 

FIG. 3. Phase diagram of a Q ID conductor in the temperature (T) vs magnetic field (H) plane 



calculated from Eq. (11) with neglected umklapp interaction (73 = 0. Various symbols denote the 



integer values of the quantized wave vector of FISDW, Qx = NG. 

FIG. 4. Temperature dependences of longitudinal resistivity normalized to the dimensionless 
umklapp scattering amplitude, pxx{T)/g3, at zero magnetic field for different values of the phase 
iph- The solid lines correspond to tj, = 0, in which case resistivity diverges at T ^ Tq, where 
To = 14.7 K is the SDW transition temperature. The dotted lines correspond to t[ = 20 K, in 
which case there is no SDW transition at zero magnetic field. 

FIG. 5. Temperature dependences of the umklapp resistivity at = 20 K for different values 
of the phase (ph in the magnetic fields 5, 15, and 25 T (curves a, b, and c, respectively) and without 
magnetic field (dots). 

FIG. 6. Temperature dependences of the umklapp resistivity at = 30 K (top panel) and 
tj, = 40 K (bottom panel) for ipf, = in the magnetic fields 5, 25, and 50 T (curves a, b, and c, 
respectively) and without magnetic field (dots). 

FIG. 7. Temperature dependences of the 2kp component of the NMR relaxation rate (in arbi- 
trary units) at zero (dots) and nonzero (solid curves) magnetic field for = 20 K. 
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